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ABSTRACT 



We examine a generalized VT -symmetric quartic anharmonic oscilla- 
tor model to determine the various physical variables perturbatively in 
powers of a small quantity e.We make use of the Bender-Dunne opera- 
tor basis elements and exploit the properties of the toatally symmetric 
operator T mn . 
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1. Introduction: 

During recent years VT -symmetric quantum mechanics has emerged as an area of high theoretical 
interest (e.g., [1-10] and references therein). For one thing, VT -symmetry is a weaker condition 
compared to the usual Hermiticity but exhibits all the essential properties of a Hermitian quantum 
Hamiltonian. For another, VT -symmetry opens up the window to the non-Hermitian world, thus 
enabling one to address a much broader class of Hamiltonians. 

Although the current interest in VT -symmetry stems from the 1998 seminal paper of Bender 
and Boettcher [1] where it was shown that for a certain class of VT -symmetric Hamiltonians the 
spectrum remained entirely real, discrete and bounded below, the concept of VT -symmetry had 
its roots in some earlier independent works as well. These include the ones of Caliceti et al [11] 
and Bessis and Zinn- Justin who studied a cubic anharmonic oscillator model with an imaginery 
coupling and that of Buslaev and Greechi [12] who analysed the spectra of certain non-Hermitian 
versions of the quartic anharmonic oscillator. 

Recently Mostafazadeh [13] has revisited the question of observables for the VT -symmetric cubic 
anharmonic oscillator problem and, in this regard, has performed a perturbative calculation of the 
physical observables including investigation of the classical limit. Motivated by Mostafazadeh's 
work, we examine, in this note, the VT -symmetric version of a generalized quartic anharmonic 
oscillator described by the Hamiltonian 

2m 2 y ' 

with (/i, v G R) that includes a cubic anharmonicity as well. 

Noting that a C -operator can be introduced [14] in the physical Hilbert space 1i v h ys subject to 
a CVT -inner product [15], and that it commutes with both H and VT , we show that for the above 
H an equivalent Hermitian Hamiltonian h can be set up. The classical Hamiltonian H c is then 
obtained in the limit h — > O.The physical position and momentum operators X and P ,which are 
actually i] + -pseudo-Hermitian for the metric operator r] + and related to the conventional position 
(x) and momentum (p) operators by the same similarity transformation that links H and h, clearly 
turns out to be VT -symmetric, a result similar to the VT -symmetric cubic oscillator. We also 
calculate the eigenvalues of H based on the first-order Rayleigh-Schrodinger perturbation theory 
upto and including terms of order e 3 . Further we determine the conserved probability density for a 
given state vector ip e 'H p h ys - It should be mentioned that our calculations are somewhat different 
from Mostafazadeh's in that we have exploited the symmetrized objects T m n [17] satisfying commu- 
tation (lowering type) and anti-commutation (raising type) relations to write down the perturbative 
expansion of the C -operator. 



2. Basic equations: 

To ensure the reality of the spectrum of a diagonalizable operator it is necessary that the 
Hamiltonian if must be Hermitian with respect to a positive definite inner product < ., . >+. The 
latter can be expressed in terms of a positive definite metric operator 77+ : 7f — > 7f of the reference 
Hilbert space H in which if acts: 

<.,.>+=< .,77+. > (2) 

where 77+ belongs to the set of all Hermitian invertiable operators 77 : 7i — > 7f satisfying iff = r/Hr]^ 1 
[18] and can be expressed as 

77+ = e-0 (3) 
where Q is Hermitian. In terms of 77+, C admits a representation 

C = V V+ = V+V (4) 

The operator C commutes with both H and VT : [C, if] — ,[C,VT] = and mimicks the charge 
conjugation operator in particle theory. 

Any Hermitian physical observable O G TL p h ys can be converted to a Hermitian operator o E Ti, 
by the transformation 

= p- l op (5) 
where p — y/r] + is a unitary operator and because of (3) may be given by 

p = e~ Q ' 2 (6) 

In view of (5) we can write 

H = p- l hp (7) 

where h is the corresponding Hermitian Hamiltonian. The classical Hamiltonian // c (x c ,p c ) is 
obtained from /i(x c , p c ) by the relation 

#c(x c , p c ) = lim ft(xc, p c ) (8) 

n— »0 

where the limit is assumed to exist. 

For the sake of convenience let us introduce a set of new variables 

X := h' 1 *, P := p, M := m 1/2 hp, e := mh 3 e (9) 

In terms of X and P, if (x, p) -> If (X, P) with 

if (X, P) = H (X, P) + eH^X, P) + £ 2 /f 2 (X, P) (10) 



H {X,P) = \P 2 + \M 2 X 2 (11) 

H 1 (X,P)=tX 3 (12) 

H 2 (X,P) = -X 4 (13) 

H(X,P) = mH(x,p) (14) 



3. Determining the Q and C-operators: 

Using (3)in (4), we consider the general form of C as 

C = e^ x ^V (15) 

It has the following properties: 

[C,VT] = (16) 
C 2 = 1 (17) 

[C,H{X,P)] = (18) 

but [C, V]^0 and [C, T] ^ 

Substitution of C from (15) into (16)implies 

e Q (X ,P) VVT = VTe Q {X ,P) p (1Q) 

showing Q(X, P) to be an even function of X: e Q(x - p) = e Q(-X ' p) .That Q(X, P) is an odd function 
of P follows from the consideration (17): 

e Q(x,p) Ve Q(x,p) v = j (20) 

which yields e Q(x ' p) = e -Q(~ x - p ) . 

We now expand of Q(X, P) in a series of odd powers of e , namely 

Q(X, P) = eQ 1 (X, P) + e 3 Q 3 (X, P) + e 5 Q 5 (X, P) + e 7 Q 7 (X, P) + 0(e 8 ) (21) 

Using(18),it follows that 

e Q(x,P) Ho _ Hoe Q(x,P) = £{e Q(x,P) Hi + Hie Q(x, P)) _ £ 2 {e Q(x,P) H2 _ # 2e Q(*,P)) 

Left multiplying both sides leads to, 



Using Baker-Campbell-Hausdorff identity, i.e, 

e~ A Be A = B + [B,A] + ±[[B,A],A] + A], A], A] + .... (23) 

we can arrange the above expression as 

-^[[[[[H ,Q],Q],Q],Q],Q]-^[[[[[[H ,Q],Q],Q],Q],Q],Q] 
-^[[[[[[[H ,Q],Q],Q],Q],Q],Q],Q] 

= 2SH, + e[H 1 , Q] + ±[[H U Q],Q] + Q],Q],Q] + M H i> Q\MMM 

+§[[[[[H 1 ,Q],Q],Q],Q],Q] + §[[[[[[H 1 ,Q],Q],Q],Q],Q],Q] 

[[[[[[[Hu Q],Q],Q],Q],Q},Q],Q}+s 2 [H 2 ,Q] + ^[[H 2 ,Q],Q] + ^[[[H 2 ,Q],Q],Q] 
+^[[[[H 2 ,Q},Q},Q],Q] + §[[[[[H 2 ,Q],Q],Q},Q],Q] 

+§[{[1{[H 2 , Q],Q],Q],Q],Q],Q] + %[[[[[[[H 2 , Q],Q],Q],Q],Q],Q],Q] (24) 
where we have taken the terms upto order e 7 

Substituting (21) into (24) and equating terms of order e,e 3 ,e 5 ,e 7 we get, 
[H ,Q 1 ] = -2H, (25) 

[HoM = -§[Qi,[Qi,tfi]] + [Qi,tf 2 ] (26) 

[HoM = -\{[QzAQi,m + [QiAQ^m) + MQ^QiAQiAQum\\ + [Q^m (27) 
[HoM = -|([Q5,[Qi,^i]] + [Q 3 ,[Q3,^i]] + [Qi,[Q5,^i]]) 

+MIQ^ [Qu IQuHM + [Qi, [Qs, [Qu [Qu h 1 ]]]] 

+[Qi, [Qu [Qs, [Qi, HM + [Q l5 [Q l5 [Q u [Qs, H 1 }}}}) 

-TsktQi, [Qi, [Qu [Qi, [Qu [QuHMW + [Qs,H 2 ] (28) 
note that terms of order e 2 , e 4 , e 6 , i.e, even powers of e give no new results. 

To solve for (25),(26),(27)and (28) we introduce, following Bender and Dunne [16] the totally 

symmetrized sum T TyS over all terms containing r-factor of P and s-factor of X. For example,we have 

T , = 1 

7\,o = P 

T 1)2 = \(PX 2 + XPX + X 2 P) 

T ,3 = X 3 

T 3 ,i = \{XP 3 + PXP 2 + P 2 XP + P 3 X) 



and so on. 



We thus get 

Qi = -lM~ 4 T 3fi -2M- 2 T h2 (28a) 

Q 3 = C-§M~ 10 - fX- 8 )T 5 , + (fM- 8 

-167W- 6 )T 3i2 + (8M- 6 - 8_M- 4 )T M - (127W- 8 - 8M- 6 )T lfi (286) 

Q 5 = (6|87W- 12 -1287W- 10 + 1287W- 8 )T 1i2 + (-64A^- 10 + 327W- 8 -327W- 6 )T 1i6 

+ (^pM- 14 - 256M' 12 + 6 -fM- w )T 3fi + (-f^AT 12 + ?fM~ w - 128M' S )T 3A 
+(- 5 -fM~ u + 128A^- 12 - 128A^- 10 )T 5i2 + (-f^AT 16 + f^M" 14 - ?fM~ l2 )T 7fi (28c) 

q 7 = (^ M - 22 -^M- 20 + ^M- 18 -^M- 16 )T 9fi 

+ (9^2^-20 _ 62^8^-18 + 34816^-16 _ l 02 4^- 14 )T 7i2 

H^WM- 18 - M- 16 + 7 -^M~ u - 1536M' 12 )T 5A 
+( J WM~ 1G - ^M' 1A + 2 -^M~ 12 - ^P_M- 10 )T 3 , 6 
+^2m M -u _ 256M- 12 + 128M- 10 - 1287W- 8 )T li8 
+ (_222|24_ M -2o + eiap^-is _ ^1 M -^ + 3584.M- 14 )T 5 , 
+(-^f§F-^" 18 + ^AT 16 - ^AT 14 + 5120Ar 12 )T 3 , 2 

+ ( _39|36_ M -16 + 40^2^-14 _ UIQM^ 12 + 1280A^- 10 )T 1)4 

+ ( 46|76_ M -18 _ 4947^-16 + 1 536 ^-14 _ 1280M~ 12 )T lfi (28d) 

Explicit forms of Q(X, P) and C(X, P) are then obtained by substituting the above expressions for 
Qi,Q 3 ,Q5,Q 7 into (21) and (15). 

4. Determining the X and P operators: 

Now we calculate the physical position and momentum operator X and P from the previously 
introduced variables X,P using the similarity transformation (5): 

X = p- 1 Xp=e%Xe~% (29) 

p = p- 1 Pp = e*Pe-% (30) 



Using (23), we obtain for X and P 

X = X + e{2iM- 4 P 2 + iM~ 2 X 2 ) + e 2 {2M~ 6 XP 2 -2tM- 6 P-M' 4 X 3 ) 

+£ 3 [(-^M- w + 16M~ 8 )iP 4 - (5M- 6 - AM~ A )iX A 

-(IfM' 8 - A8M' 6 )XP + (f AT 8 - 24M~ 6 )iX 2 P 2 

+{fM- 8 - lQM- & )i] + 0{e A ) (31) 

P = P - e(2iM- 2 [XP - D) + e 2 (2M- 6 P 3 - M~\X 2 P -iX)) 

-i£ 3 [(16M- 8 - 16M- 6 ){XP 3 - \iP 2 ) 

+{\QM- & - \QM- A )(X 3 P - \iX 2 )} + 0{e A ) (32) 

where we retained terms of order of e 3 . 
From (31) and (32), we easily see that 

VX.V ^ -X 
VW ^ -P 
TXT^X 
TPT ^ -P 

but 

VTX.VT = -X 

VTPVT = P (33) 

From (33) we conclude that the physical position and momentum operator, i.e, X and P G H p h ys 
are PT-symmetric. 

' 5. The equivalent Hermitian Hamiltonian: 

For the operator p, the corresponding Hermitian Hamiltonian h(X, P) is given according to 
(7)with previously introduced variables X,P as 

h(X, P) = e- Q/2 H{X, P)e Q ' 2 (34) 
Introducing the perturbative expansion for h(X, P) as 

oo 

h = h^e* (35) 

i=0 



and using (10),(21),(25)-(28) and (34) along with the Baker-Campbell-Hausdorff identity(23), we 
obtain for the various coefficients ,i=0, 1,2,3,...., the results 
= H 

h& = H 2 + l[H 1 ,Q 1 ] 

= 1[H 1 ,Q 5 ]-^([[[H 1 ,Q 1 },Q 1 IQ 3 } + [[[H 1 ,Q 1 },Q 3 IQ 1 ] 

+[[[H 1 ,Q 3 ],Q 1 ],Q 1 }) + ^[[[[[H u Qi],Qi],Qi],Qi],Qi] (36) 
with the odd ones vanishing: = ,i=l,3,5,7, 

Keeping terms up to the order e 5 ; h(X, P) can thus be expressed as 

h(X,P) = |(T 2i0 + ^ 2 T 0i2 )+£ 2 [-|^- 4 + (|^- 2 -l)r ,4 + 3^- 4 r 2i2 ] 

+e^[-(36M- w - 247W- 8 )T 4 , 2 + (27M- 10 - 24M- 8 )T 2fi - (§ M~ 8 - 36M- 6 )T 2A 
+(^M- 8 - 12^- 6 )T 0i2 - (fM- 6 - 6^- 4 )T 0i6 + 2M- 12 T 6fi ] + 0( £ 6 ) (37) 

If now we consider the normalized eigen vector \n > of the conventional Harmonic oscillator Ho, 
then we can easily calculate E n for H by the first order Rayleigh-Schrodinger perturbation theory. 
We obtain upto the terms of order e 3 

E n = M{n+\)+e 2 <n\h^\n> +0(e 4 ) 

= M{n + \) + ^[^(30n 2 + 30n + 11) - (6n 2 + 6n + 3)] + 0(e 4 ) (38) 

6. Classical Hamiltonian: 

Employing (9) , we have, 

T r , s = h- s S r , s (39) 

where S r>s be the totally symmetrized sum of all terms containing r-factor of p and s-factor of x. 
Specifically 

<So,o = 1 
So,i = x 
S 3 ,o = P 3 

Si,3 = \{x 3 p + xpx 2 + x 2 px + px 3 ) 

and so on. 

Now from (14), 



Keeping terms up to the terms of order of e 3 in (37) and using the relations (9), (39), (40) we finally 
obtain 

h(x,p) = ^ + -^V + e 2 m[(-m-V~ 2 -/i 2 )^ 4 -2m- 2 ^V 4 -6^- 2 ^-V-3m-V"W] + 0(e 4 ) 

2 2 

(41) 



The corresponding classical Hamiltonian can be read off from (8): 
H c (x c ,p c ) = lim ft ^o h(x c ,p c ) 

where 

MM = T^Fw c (43) 

A position-dependent mass M(x c ) is implied by (43) for the classical particle whose dynamics is 
dictated by the Hamiltonian H c (x c ,p c ). 

7. Conserved probability density: 

For a given state vector -0 e TC p h ys , the perturbation expansion for the corresponding physical 
wave function is 

ty(x) = < x | e - ^ | -0 > 



_ <T iv» i > 



13. 



= ^(x)+ < x I -'f I > e+ < x | ^ | if; > e 2 + < x | (-Sk - ^) | ^ > £ 3 + 0(e 4 ) (44) 

Using (9), (28a) and (28b) we obtain from (44) 

= (1 + eLi + e 2 L 2 + e 3 L 3 )^(x) + 0(e 4 ) (45) 

where 

Li = -f^Qi (46) 
^2 = f % % Qi (47) 
L 3 = -r£ft 9 Qs - ^ft 9 & 3 (48) 



48 

and, 



Qi = -|m- 2 ^-V" 4 ^3,o - 2m- 1 h- 2 n- 2 S h2 (49) 
Q 3 = (f*m- 5 7rV" 10 -f™- 4 ft~V 8 )S 5 ,o + (^^ 



Employing (46)-(50) into (45), we find the conserved probability density g associated with a 
given state vector ip e H-phys as 

e (x) = N' 1 | | 2 (51) 

where 

/oo 
I | 2 (52) 
-oo 

8. Conclusion: 

We have carried out a perturbative treatment to study a -symmetric quartic anharmonic 
oscillator model. We have shown possible to set up an equivalent Hermitian Hamiltonian by em- 
ploying a similarity transformation. Such a Hamiltonian has a classical limit too. Physical position 
and momentum operators have been determined perturbatively and energy eigenvalues are obtained 
in the framework of first order Rayleigh-Schrodinger perturbation theory. In all these calculations 
we have kept terms up to and including those of order e 3 . The conserved probability density is 
also determined. Finallylet us mention that in the absence of the quartic term in (1) our results 
essentially reduce to those of Mostafazadeh's [13] for the physical observables X and P, equivalent 
Hermitian Hamiltonian h(X, P) and the energy spectrum derived from the first order Rayleigh- 
Schrodinger perturbation theory. 
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